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Abstract We review recent theoretical models and ex- 
periments dealing with the diffraction of neutral atoms 
by a reflection grating, formed by a standing evanes- 
cent wave. We analyze diffraction mechanisms proposed 
for normal and grazing incidence, point out their scopes 
and confront the theory with experiment. 
PACS: 32.80.Lg Mechanical effects of light on atoms - 
42. 25. Fx Diffraction and scattering - 03.75.Dg Atom in- 
terferometry 



' 1 Introduction 

[ An evanescent wave is the light field formed in vacuum 
• above a dielectric when a light beam undergoes total in- 
' ternal reflection at the vacuum-dielectric interface. The 
evanescent wave propagates parallel to the interface and 
decreases exponentially with distance from the interface, 
at a scale sligthly smaller than the optical wavelength. 
In 1982, Cook and Hill |§ proposed using this spatially 
inhomogeneous light field to construct a repulsive optical 
potential able to reflect neutral atoms with normal ve- 
locities (perpendicular to the interface) of several cm/s. 
In 1989, Hajnal and Opat |^ proposed to combine two 
counterpropagating evanescent waves, creating a mirror 
with a spatially periodic modulation in order to realize 
a reflection grating (cf. figure 0). At grazing incidence, 
the diffraction angles are greatly enhanced similar to X 
ray diffraction from optical gratings, and one obtains a 
beamsplitter useful for atom optics ||^,^ and interferom- 
etry |. 

It took several years of theoretical and experimen- 
tal work 1^,^ before efficient diffraction from a station- 
ary evanescent wave was observed This is due, 
on the one hand, to the experimental difficulties: the 
diffraction angles are only in the mrad regime because 



of the small atomic wavelength compared to the grating 
period; the grating has to achieve both reflection and 
diffraction of the atomic beam; at grazing incidence, the 
coupling between the diffraction orders is reduced due to 
the Doppler effect. Diffraction was first achieved exper- 
imentally when all these were overcome: in 1993, in the 
experiment at Bonn university (Germany), the atomic 
beam had been slowed down while in 1996, the group 
at Paris-Nord university (Villetaneuse, France) tilted the 
beam with respect to the evanescent wave's propaga- 
tion direction At normal incidence, the ENS group 
(Paris, France) reported diffraction in the time domain 
(1994) fill, and in 1996, the group at the Orsay Institut 
d'Optique (France) observed that the evanescent wave 
yields a nonspecular (diffuse) reflection , washing out 
the diffraction peaks, unless in 1997 dielectric surfaces 
with high surface quality were used 0. 

On the other hand, the theory of the atomic reflec- 
tion grating encountered a number of difficulties: for in- 
stance, the atomic motion is very different from a trans- 
mission grating because the atoms are also reflected from 
the grating, excluding the elimination of one direction of 
motion in the Raman-Nath approximation. In 1992, the 
group at Bonn university solved the Schrodinger equa- 
tion for a two-level atom in a coupled-wave formalism 
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Fig. 1 Schematical setup of an evanescent wave reflection 
grating, (a): normal incidence, (b): grazing incidence. 
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p3t and showed that atomic diffraction is possible if 
the Doppler shift is comparable to both the optical po- 
tential and the detuning. The proposed diffraction pro- 
cess involves a delicate interplay of adiabatically fol- 
lowed dressed states and nonadiabatic transitions be- 
tween them. This approach has been used to interpret 
numerical calculations and allowed a quanti- 

tative understanding of experiments on Doppleron reso- 
nances 0. 

For the case of normal incidence and large detun- 
ings, the group at Institut d'Optique (Orsay, France) 
developed in 1994 an alternative theory based on scalar 
diffraction, similar to light diffraction from thin phase 
objects This approach showed that the atomic 

wave is very sensitive to spatial modulations of the 
evanescent wave because of its small wavelength, in 
agreement with the Orsay experiments at normal inci- 
dence [|o[|l^. At grazing incidence however, the same 
theory predicts vanishing diffraction, in contrast to ex- 
periment and the two-level theory jl^ . 

In 1995 the Canberra theory group realized that 
a numerical integration of the coupled wave equations 
for a two-level atom at grazing incidence does not repro- 
duce the Bonn experiment |^ , indicating that Doppleron 
resonances do not provide the relevant diffraction mech- 
anism. Since 1996, Gordon and Savage Henkel et al. 
pO| and Deutschmann |^ pointed out that the observed 
diffraction at grazing incidence can only be explained in 
terms of Raman transitions between Zeeman sublevels 
of the atomic ground state. A similar mechanism had 
been studied in 1993 by the Bonn group |2|] in view of 
building a grazing incidence beamsplitter from a running 
evanescent wave and a static magnetic field. 

The present paper aims to develop a unified picture 
of these different theoretical approaches: we review the 
diffraction mechanisms mentioned above and compare 
them to each other and to the diffraction experiments 
reported so far We start in Sec^jfrom the most 

simple description in terms of a 'one-level' optical poten- 
tial and discuss the results of a perturbative calculation 
for the diffraction pattern. For normal incidence diffrac- 
tion, this theory has been extended using semiclassical 
techniques to compute higher-order diffraction peaks. 
We then pay particular attention to the case of grazing 
incidence. In Sec.^, the two-level model of Deutschmann 
et al. is reviewed. We compare it to the one-level model 
and give an interpretation why these models make dif- 
ferent predictions for the diffraction probabilities. We 
then apply the two-level model to the diffraction exper- 
iments and give an analytical estimate of the Doppleron 
coupling, showing it to be too small to account for the 
observed diffraction patterns. On the other hand, the 
situation is quite different for Raman couplings, as we 
discuss finally in Sec.^. A simple condition for the opti- 
mum diffraction efficiency is derived and we outline some 
perspectives for experiments with spin-polarized atoms. 



2 Diffraction of one-level atoms 

2. 1 Optical potential 

In coherent atom optics, one frequently works in con- 
ditions where the atom is driven non resonantly and 
at low saturation, in order to avoid spontaneous emis- 
sion processes. Under these circumstances, the theoret- 
ical description of the atom-laser-interaction simplifies 
because the excited state may be eliminated adiabati- 
cally. The atomic dynamics is governed by a Schrodinger 
equation for the ground state wavefunction only ('one- 
level-atom'), where the light field enters via the so-called 
'optical potential' or dipole potential 

T/opt(r) = ^|E(r)p (1) 

In this expression, d is the reduced matrix element of the 
atomic dipole operator, 5 = lol — is the laser detun- 
ing, lol {^a) is the laser (atomic resonance) frequency, 
respectively, and E(r)e~*'^^* -I- c.c. is the laser field. 

To begin with, we consider only a single state in the 
atomic ground state. The atomic wavefunction ^{v, t) 
and the Schrodinger equation then become scalar quan- 
tities (without Zeeman sublevel indices) 

'^^'"^ = -^^^^ + ^opt(r)V (2) 

The AC-Stark-shift has now acquired the meaning of a 
ponderomotive potential for the center-of-mass motion 
of atoms without internal structure. Such an approach 
has found widespread use in theoretical investigations 
of scalar atom optics, e.g. for the study of transmis- 
sion gratings (see [|] for further references). Although 
this approach is historically not the first for the reflec- 
tion grating [the first investigations (|,|l^ were done for 
a two-level atom at arbitrary saturation, see Sec.^, it 
gives the most intuitive understanding of the problem. 
Its most straightforward application is the description of 
'thin' optical elements, i.e. elements which are character- 
ized by the dynamical phase the atoms accumulate when 
passing the optical element. Moreover, it allows one to 
develop far-reaching analytical methods and results that 
prove useful in the interpretation of more complicated, 
even non-scalar models. 

2.2 Kinematics 

2.2.1 Reflection grating. In contrast to a usual trans- 
mission grating the atomic reflection grating has to 
achieve two functions, diffraction and reflection, by a 
suitably shaped optical potential. A reflective optical 
potential with a sharp intensity variation along one di- 
rection and a periodic intensity modulation along the 
perpendicular direction, that does this job, can be real- 
ized by overlapping two evanescent waves with wavevec- 
tors Qe^ and —Qe^ parallel to the dielectric surface 
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Fig. 2 Optical potential (g) of a partially stationary evanes- 
cent wave (contrast e = 0.3). 



{Q > k = ujl/c). Denoting the corresponding field am- 
plitudes by E±(r) = E± exp (ztzQa; — the optical 
potential (|l|) becomes 



with 



Vopt{x, z) = Vmax(l + ecos2Qx)e~ 



Kna.= ^(|E+P + |E_n 

2ReEl • E_ 



|E+|2 + |E_ 



(3) 

(4) 
(5) 



The overall height of the potential is given by Vniax, 
while the dimensionless quantity e determines the con- 
trast of the stationary wave. The potential (|^) is plot- 
ted in FigH The period of the optical potential is equal 
to tt/Q = A/(2nsin0i) and typically a fraction of the 
optical wavelength (n: index of refraction of the dielec- 
tric). The same is true for the decay length 1/2k — 



A/[47r( 



n sm 



_ 1)1/2] 



Atoms impinging on this optical potential will experi- 
ence a reflection of the z-component of their momentum 
if their kinetic energy in this direction does not exceed 
the potential height and they can pick up momentum 
changes ±2Q from the periodic modulation. 

In Eq. , we have neglected the (typically gaussian) 
intensity profile of the evanescent wave, assuming plane 
waves incident inside the dielectric. The present picture 
also neglects the van der Waals interaction between the 
atom and the dielectric This interaction is attrac- 
tive and becomes larger than the optical potential for 
distances smaller than roughly \/Ati. The full potential 
is well approximated by the optical potential (|^) for dis- 
tances larger than this value. 

2.2.2 Diffraction channels. We now introduce the 
diffraction channels for an atom incident from infinity 
onto the stationary evanescent wave. From energy and 
momentum conservation for the asymptotic states (at 
large distance z — > oo from the grating), one finds that 



the diffracted matter waves are characterized by the fol- 
lowing wavevectors 



h 



2nQ, n = 0,±l,±2,. 



kzn = ^Jkli - AnQ{kxi + nQ) 



(6) 
(7) 

(8) 



where the incident wavevector equals {k^i, kyi, —kzi). In 
the direction parallel to the grating (the x-direction), 
the atomic wavevector changes by an integer multiple of 
the 'grating vector' 2Q. By translation invariance, the y 
component is conserved. Finally, the normal (or z-) com- 
ponent is obtained from energy conservation, the opti- 
cal potential being time-independent. These wavevector 
transfers may also be understood in the photon picture: 
in the diffraction process, the atom absorbs a photon 
from one travelling evanescent wave and emits a (stim- 
ulated) photon into the other wave, thereby returning 
to the ground state. This can be repeated n times for 
the nth diffraction order. The momentum exchange dur- 
ing the process leads to Eqs.(^, |^). Energy is conserved 
[Eq.(H)], corresponding to an elastic scattering process, 
because the two light waves have the same photon energy 
(frequency) . 

From Eqs.(||-H), the diffraction angles (with respect 
to the surface normal) are easily obtained through 
tan0„ = kxn/kzn- In this article, we focus throughout 
on two limiting cases: normal and grazing incidence. At 
normal incidence and under typical experimental condi- 
tions, the first term in the square root in Eq.(|^) dom- 
inates the two others: the atoms being incident with a 
velocity much larger than the recoil velocity hk/M, their 
incident wavevector k^i is much larger than the optical 
wavevector k Q. As a consequence, the normal com- 
ponent of the atomic wave vector changes negligibly, and 
the diffraction angles are small (of order 6'„ ~ 2nQ/kzi). 
The diffracted beams have a transverse velocity of the 
order of a few recoil velocities. 

At grazing incidence, the second term in the square 
root, —AnQkr^i^ is comparable to the first, and diffraction 
leads to a substantial change of the normal wave vector 
component. The diffraction angles are hence compara- 
ble to the angle of incidence, and the diffraction orders 
are conveniently separated in space even for a quite fast 
atomic beam. This enormous enhancement of the diffrac- 
tion angles represents the main motivation for atomic 
diffraction at grazing incidence. 

Note that the above effect may be interpreted in an 
alternative way in the reference frame co-moving with 
the atomic beam. In this frame, the frequencies of the 
two evanescent waves differ by twice the Doppler shift 

= Qvxi — ^Qkxi/M, and diffraction is accompanied 
by an energy transfer 2nh5u- With this reasoning, one 
also finds Eq.(^), up to the 'recoil shift' term An^Q^ in 
the root (usually a small correction compared to the first 
and second terms). 
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2.3 Diffraction intensities: Born approximation 

2.3.1 Diffraction amplitudes. Up to now, we have only 
specified the kinematics of the diffraction process. In 
this subsection, we present a first approach to calcu- 
late the diffraction pattern^ i.e. the fraction of atoms 
diffracted into the different orders. To this end, we in- 
troduce the asympotic form of the atomic wavefunction 
in the region of vanishing potential (a time-dependent 
factor exp{~iEt/h) has been separated) 

z -> +0O : (9) 
ip{x, z) = expi(/c„x - kziz) -I- ^ a„ expi(fca,„a; -I- k^nz) 

n 

The first term represents the incident wave and the sum 
the diffracted waves. The coefficients a„ are the diffrac- 
tion amplitudes that are related to the diffraction prob- 
abilities 7r„ according to 



.Refer 



(10) 



[Diffraction orders with imaginary kzn are evanescent. 
They are localized in the optical potential and do not 
contribute to the wavefunction at infinity.] 

To calculate the diffraction amplitudes a„, one has, 
in principle, to solve the stationary Schrodinger equa- 
tion subject to the asymptotic condition (|^), expand the 
wavefunction in the region z —> +oo into a Fourier se- 
ries with respect to x and read off the expansion coeffi- 
cients. This is a difficult task that may only be accom- 
plished numerically, using e.g. a coupled-wave expansion 
(see Sec-IJand Refs.[|3[|l|,|l|,|l§). An additional com- 
plication arises from the fact that the atoms are not 
only diffracted, but also reflected from the optical poten- 
tial. This translates into a second asymptotic condition, 
namely that the wavefunction matches an exponentially 
decreasing solution inside the potential. [Corrections to 
this picture due to tunneling through the optical poten- 
tial towards the dielectric are discussed in |2^.] For the 
same reason, we cannot use Fermi's Golden Rule to com- 
pute the diffraction amplitudes because in this formula, 
initial and final states are approximated by the corre- 
sponding asymptotic states (the plane waves of Eq.(^)) 
that are, however, a very bad approximation to the true 
wavefunction deep inside the potential. 

2.3.2 Distorted-wave Born approximation. The way to 
circumvent this problem and to obtain analytical results 
for the diffraction amplitudes dates back to the 1930's 
[^ when Lennard-Jones and Devonshire calculated the 
diffraction of a thermal atomic beam from a periodic 
(crystalline) surface. The idea is to treat first the unmod- 
ulated potential (that provides the momentum transfer 
for the basic reflection process) in an exact manner, and 
then to use perturbation theory for the modulated part 
of the potential. One obtains a modified Fermi's Golden 



Rule where wavef unctions 'tjjij(z) are used that decrease 
inside the potential and asymptotically contain an inci- 
dent and a specularly reflected wave. The matrix element 
evaluated in the Golden Rule thus reads 



+ OC 



Af 



= (V/(z)| e-^«- = / dzrf{z)e-'^'Uz) 



(11) 

This approach to the diffraction problem is called the 
'distorted wave Born approximation' (DWBA) in the fol- 
lowing. It is a perturbative one and of course valid only 
for a stationary evanescent wave with low contrast such 
that nonzero diffraction orders are weakly populated. 

When the DWBA is applied to the optical poten- 
tial (||) , one notes first that the Schrodinger equation for 
the flat potential (an exponentially increasing potential 
barrier) has an analytical solution j2^,^, and second 
that the matrix element of the modulated part of the po- 
tential may be evaluated in closed form [|2^,^. In the 
present notation, one thus obtains the following result 
for the diffraction probabilites [neglecting exponentially 
small corrections of order exp(— Trfczi/re)] 

2 

P^[{k,^±i-k,,)/K\ 



7''±l,DWBA 



+ k. 



2k 



(12) 



(13) 



where the dimensionless factor /3(^) is defined by 

sinh(7r^/2) 

and reaches its maximum value (unity) for f = 0. 

Note that the DWBA gives the lowest order contri- 
bution to the diffraction pattern. It hence determines 
only the intensities of the first-order diffraction peaks 
{n = ±1). Calculations to second order are possible, 
though much more involved because one has to integrate 
over a continuous spectrum of intermediate states. 

2.3.3 Physical interpretation. 

Diffraction efficiency. We first note that the diffraction 
process may be quite efficient even for a small contrast 
e. Indeed, if the factor /? in ( p^ ) is close to its maximum 
value, the diffraction peaks have a height proportional 
to the square of the product efczi/t = 27re cos ^^/(KAdB), 
and since AdB, the wavelength of the incident matter 
wave, is typically much smaller than the decay length 
1/k, the contrast e is multiplied by a large number. 

To understand this feature, let us consider the po- 
tential of Figj^. Obviously the contours of the poten- 
tial equalling a given atomic kinetic energy are nearly 
sinusoidal curves in the xz plane. This suggests to re- 
place the optical potential of the stationary evanescent 
wave by an infinitely high wall with a spatially modu- 
lated position [the 'corrugated hard wall' potential used 
in atom-surface scattering theory |30t|3l|| ] 

(14) 



z{x) — (e/2/v) cos2Qx. 
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The atomic wave acquires a position-dependent phase 
shift S4>(x) = 2kziz(x) upon reflection from this bar- 
rier. The phase of the outgoing wave hence shows a 
phase modulation u cos 2Qx with a modulation depth 
u = ekzi/n- The phase-modulated wave may be ex- 
panded into sidebands, and the intensity of the first- 
order sidebands is given by (for small contrast) 



hard wall: 



- ) (15) 



This coincides with Eq.(|l2|) in the regime kz,±i ~ kzi 
where (3=1. The intensity of the sidebands is compa- 
rable to that of the carrier when the phase modulation 
depth approaches unity. This means that the amplitude 
of the hard-wall corrugation, e/2K, must be comparable 
to the wavelength of the incident wave. The latter being 
much smaller than the optical wavelength, a low contrast 
e is sufficient to deeply modulate the atomic phase. 

This picture of course ignores the finite width of the 
optical potential. We show now that this width is related 
to the behaviour of the factor f3 in ( |l2|) as a function of 
the normal wavevector transfer Ak^n — k^n ~ kzi- 

Normal vs. grazing incidence. The function (3{Akznl n) 
( p^ ) rapidly vanishes as soon as the normal wavevector 
transfer Akzn exceeds the evanescent wave decay con- 
stant AC. It is close to its maximum value of 1 near normal 
incidence where the normal wavevector transfer Akzn is 
small and varies quadratically with the grating vector Q 
[Eq.(||), recall that typically Q <C kzi]. On the contrary, 
at grazing incidence, Akz,±i is much larger than k and 
(3 decreases exponentially fast: 



|Z\fc,|>K: (3{Akz) 



n\AK_ 



■ exp 



TT\Akz 

2k 



(16) 



This fact practically induces a cutoff in the diffraction 
process for momentum transfers \Akzn \ ^ k. In particu- 
lar, the DWBA predicts a vanishing diffraction efficiency 
for one-level atoms at grazing incidence. This feature is 
in contradiction with the experiments at grazing inci- 
dence [^|9|, and we would like to discuss it in more de- 
tail. 

Consider the matrix element Afi ( [III ) that appears 
in Fermi's Golden Rule when the DWBA is used. The 
z-integral is in fact limited to a narrow interval of width 
1/k because in the classically forbidden region z — > — cx), 
the wavefunctions ipij{z) vanish exponentially whereas 
in the region z — > -f oo, the exponential potential limits 
the integral. The integrand in Eq.(pl]) is hence an oscil- 
lating function due to the interference between the initial 
and final waves, with an envelope of width ^ Ap- 
plying an argument familiar from the theory of Fourier 
transforms, we may now estimate that Afi is signifi- 
cantly different from zero only if the difference between 
the wavevectors kzi^ kzf is smaller than the inverse width 
of the envelope: 



Otherwise stated, the diffraction from a standing evanes- 
cent wave only provides a 'photon momentum' of order 
TiK in the normal direction.^ (For an atomic transmis- 
sion grating, this was pointed out in 1987 by Martin et 
al. 1 32 .) The functional form Afi depends on the enve- 



lope of the integrand in the matrix element. Noting that 
this envelope is a smooth function of position, the gen- 
eral properties of Fourier transforms imply that A fi be- 
comes exponentially small for wavevector transfers much 
larger than the limit (|l7| ); this may also be checked from 
the exact result ([l^ , |l6| ). This behaviour is in sharp con- 
trast to the hard-wall potential that efficiently provides 
large wavevector transfers in the normal direction, as 
seen in the model of the preceding paragraph. In their 
paper |2^, Armand and Manson used the modulated 
exponential potential in order to study the impact of a 
finite-width interaction potential in atom-surface scat- 
tering, in comparison to the corrugated hard wall model. 
They also observed that diffraction decreases if the po- 
tential gets 'softer'. 

We are hence led to the conclusion that in the frame- 
work of the present (scalar) model, no atomic diffraction 
is predicted at grazing incidence, which indicates the 
need for an alternative approach in order to interpret 
the experimental results [0,^ obtained in this geometry. 
This will be done in Secs]| and ^. 



2.^ Diffraction intensities: semiclassical perturbation 
method 

In this paragraph, we continue discussing the case of nor- 
mal incidence in order to compute higher-order diffrac- 
tion peaks. In fact, the sensitivity of diffraction to weakly 
modulated potentials calls for a generalization of the ap- 
proach discussed so far, to go beyond the first-order cal- 
culation of the DWBA. The method we present here 
was developed by Henkcl et al. ||l^ , elaborating on ideas 
of Cohen- Tannoudji ||l^. It still leads to analytical re- 
sults for the peak intensities and provides a simple physi- 
cal picture for the diffraction process, substantiating the 
heuristic approach of the corrugated hard-wall potential 
introduced above. The main difference to the DWBA is 
that the new approximation scheme remains valid for 
much larger values of the contrast e where the diffrac- 
tion pattern typically contains a large number of peaks. 
Another important difference to the DWBA is to use 
right from the start semiclassical concepts like trajec- 
tories and phases to compute the diffraction pattern. 
This approach is justified by the experimental condi- 
tions where the atomic wavelengths are generally much 



Af, / ^ 



< K 



(17) 



This is not in contradiction with the fact that the reflec- 
tion reverses the incident momentum, since reflection is due 
to the flat potential that is treated exactly (to arbitrary or- 
der). Diffraction, on the contrary, is caused by the modulated 
part of the potential that, to lowest order, has to provide the 
momentum transfer between the diffraction channels. 
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smaller than the optical wavelength (the typical scale of 
the diffraction grating).^ 

2.4-1 Outline. It is instructive to recall the example of 
light diffraction from an acoustic wave with a weak in- 
dex modulation, as discussed in Ch. 12 of Born and Wolf 
pSf . After the passage through the acoustic wave, the 
light field has been phase-modulated, and the diffrac- 
tion intensities are obtained from a Fourier transform of 
the field amplitude at the exit of the interaction zone. 
A simple way to calculate the optical phase modula- 
tion is to accumulate the refraction index along straight 
rays through the acoustic wave. Note that this 'recipe' is 
less accurate than the geometrical optics approximation 
since it discards the deflection of the rays due to the in- 
dex gradient. If this deflection is small, the grating may 
be called 'thin' and for this reason, the outlined method 
has been termed the 'thin phase grating approximation' 
(TPGA). Its results are identical to those obtained in 
the Raman-Nath-approximation [^,^. 

At first sight, the TPGA does not seem to apply 
to the atomic reflection grating we are interested in, 
since the 'rays' (classical paths) are substantially dis- 
torted even for the simple reflection. We now show that 
the TPGA may nevertheless be used (see [0|l^ for de- 
tails). The basic idea is to compute, for a weakly mod- 
ulated diffraction grating, the phase shift of the atomic 
wave to lowest order in the contrast e. [This approxima- 
tion does not necessarily imply only first-order diffrac- 
tion because the diffraction pattern will contain several 
orders if the resulting phase modulation depth is large.] 
In this calculation, we may use some basic concepts of 
Lagrangian classical mechanics. Recall that for a given 
trajectory linking two points, the atomic phase is equal 
to the classical action integral Sd, divided by Planck's 
constant h. According to the principle of least action, the 
classical trajectory that solves the equations of motion 
corresponds to an extremum (actually a minimum) of 
the action integral. This property allows us to compute 
the action Sd without explicitly solving the equations of 
motion: we obtain a value for Sd accurate to first order 
even if the action integral is calculated along a nearby 
path. For the evanescent wave diffraction grating with 
a weakly modulated potential, we thus get the atomic 
phase modulation from an integral along the classical 
path reflected at the flat (non modulated) optical poten- 
tial because this path is close to the 'real' path. Using 
the TPGA in this way, one obtains a simple analytical 
result because both the classical trajectory in the flat 
potential and the relevant integral along it are explicitly 

^ This remark also shows why the DWBA is useful in its 
own right; namely, it also allows one to study diffraction in 
the very-low-energy domain where the atomic wavelength 
is comparable to the scales of the diffraction grating. This 
regime has not yet been thoroughly explored in experiments, 
though it should lead to interesting quantum effects ]33|,|34[. 
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Fig. 3 (courtesy of A. Landragin, Ref.[p^) Atomic momen- 
tum distributions after normal-incidence reflection from an 
evanescent wave with a weak stationary component (inten- 
sity ratio of order 10~*). The amplitude of the standing wave 
increases from (a) to (e). The solid line is a theory based upon 
the thin phase grating approximation: squared Bessel func- 
tions convoluted with the specular reflection pattern (a) , the 
modulation index being the only flt parameter. 

known [^^,^, while the real path is more difhcult to 
calculate. 

The perturbative approach just outlined is called the 
'thin phase grating approximation' (TPGA) by anal- 
ogy to the method used in light diffraction. The atomic 
diffraction grating is 'thin' if the actual classical paths 
are only slightly perturbed by the modulated potential. 
In particular, one has to avoid focal points and caustics 
inside the interaction region (see |Q for a detailed dis- 
cussion). We note that in neutral atom scattering from 
crystalline surfaces, the TPGA is well known under the 
name 'trajectory approximation' |^,^. In the field of 
cold neutron wave optics, similar approaches have been 
discussed by Felber et al. [|8|. 

2.4-2 Discussion. When the thin phase grating approx- 
imation is applied to the atomic diffraction problem, the 
reflected atomic wave is phase-modulated with an am- 
plitude 

WTPGA = e— /3[(2g/K) tan 0,] (18) 

K 

where the cutoff function /S defined in Eq. (p^) has been 
used. The diffraction pattern predicted from the TPGA 
is thus familiar from phase modulation: the diffraction 
peaks have weights given by the squares J^iu) of Bessel 
functions (see also Fig.H). In the perturbative regime 
li <C 1, the pattern consists of a strong 'carrier' (zeroth 
order peak) plus symmetrical sidebands (diffraction or- 
ders n = ±1) with a weight of order In the opposite 
regime u ^ 1, sidebands up to order n « ±u are pop- 
ulated, with a broad maximum at the extreme orders. 
For a modulation index m w 2, the orders rt = ±1 are 
maximized to a height « 35 % each. 
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In order to establish the vahdity of the TPGA, let 
us consider the perturbative regime because there, the 
DWBA provides us with a 'benchmark result'. At nor- 
mal incidence, = and we recover, for a weak phase 
modulation u < 1, the result of the DWBA pro- 
vided the normal wavevector difference Akz^±i is small 
compared to k. Expanding this difference to second or- 
der in Q, we hence find a condition of validity of the thin 
phase grating method^ 



\Ak 



z,±l\ 



2Q2 



< K, 



(19) 



a condition that is typically satisfied in diffraction ex- 
periments. The TPGA hence extends the DWBA result 
to the regime of large atomic phase modulation. Its pre- 
dictions agree with diffraction experiments at normal in- 
cidence, as illustrated by Fig.|| taken from Rcf . [p^ . 

We conclude this paragraph with an alternative in- 
terpretation for the cutoff of the diffraction efficiency at 
grazing incidence, based on an argument given in 1987 
by Martin et al. for a transmission grating In the 
TPGA, one gets the atomic phase modulation by ac- 
cumulating the sinusoidal part of the potential during 
the reflection from the flat potential. At grazing inci- 
dence, it happens that the classical path passes through 
a large number of standing wave periods during the in- 
teraction time T. As a consequence, the standing wave 
exerts a rapidly oscillating potential that tends to aver- 
age out. Since this oscillation frequency is given by twice 
the Doppler shift 2Qvxi^ we may expect the modulation 
index to be very small in the regime 2QvxiT ^ 1. Tak- 
ing an interaction time r = l/(KWzi)) for the evanescent 
wave reflection grating, we recover the argument of the 
cutoff function [3 in Eq. (|l^) , and the asymptotic expan- 
sion ( [l6| ) confirms our expectation of a vanishing mod- 
ulation index. This cutoff has been observed in experi- 
ments performed in the time domain see also |p9| , 

f|, and very recently in spatial diffraction experiments 
y . Balykin et al. proposed a method to restore efhcient 
atomic diffraction at grazing incidence Q : they reduce 
the grating vector 2Q by using a modulated evanescent 
wave above a structured surface with a period in the /J,m 
range. 



3 Diffraction of two-level atoms 

The results of the preceding Section suggest that the 
one-level model cannot explain evanescent wave diffrac- 

^ Note that the TPGA method does not correctly de- 
scribe the wavevector transfer normal to the grating. Roughly 
speaking, this failure is related to the choice of the classi- 
cal trajectory along which the phase shift is computed. In 
the context of atomic interferometry, Borde has presented a 
generalization of the TPGA that correctly includes the term 
quadratic in Q (the recoil term) [Q. Similar improvements 
have also been reported for thermal atom scattering from 
crystalline surfaces (see |^ for a review). 



tion at grazing incidence because the 'bandwidth' of this 
grating is interaction-time limited. In order to under- 
stand why experiments at grazing incidence neverthe- 
less showed diffraction, more complex models are needed 
and have been introduced, in fact, prior to the one-level 
model discussed so far. These models have in common 
that they include, to a certain extent, the atomic mul- 
tilevel structure. Historically the first approach to do 
so is the extensively used two- level atom P,[l3|, ^[r^ , 
|6|,|4^ that allows an accurate description of the atom- 
light interaction at any degree of saturation. The 'next 
generation' also includes the magnetic sublevel struc- 
ture of the atomic ground state |l^,^,|2|. We discuss 
in this Section the two-level diffraction theory developed 
by Deutschmann et al. ||l3|. Its relation to the one- level 
model is analyzed in some detail, since it takes a quite 
different route to explain the diffraction process. 



3.1 General 

We focus again on the coherent interaction between a 
two-level atom and a monochromatic laser field and as- 
sume again that spontaneous emission is negligible. The 
atom is then described by a two-component wavefunc- 
tion (?Ag, V'e)^ whose evolution is governed by the follow- 
ing (r.w.a.) Hamiltonian matrix 



H 



2M 



f dE*{r) 
\dE{r) hS 



(20) 



The notations are identical to Eq.(|^), and we recall the 
electric field amplitude 

E{r) = {E+cjcpiQx + E-exp{-iQx))e-''''. (21) 

The laser field couples the ground and excited states, 
and therefore new diffraction channels open up where 
the atom leaves the interaction region in the excited 
state. From momentum conservation, these channels cor- 
respond to an odd number of exchanged photon mo- 
menta along the x-axis: k^i, — k^i + vQ with = 2n+ 1. 
For the odd diffraction orders, the detuning 5 between 
the photon energy and the atomic internal energy has to 
be taken into account since the atoms end up in a dif- 
ferent internal state (the scattering is inelastic). Using 
energy conservation, the normal wavevector component 
is obtained as 



odd: 



(22) 



kzi, ^ Jkl^ + 2M5/fi - t/Q(2fc„ + vQ) 



For even v, k^v is still given by Eq.(|h with n — v 12. 
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12 3 4 

distance (kz) 

Fig. 4 Adiabatic potentials in a stationary evanescent wave. 
The potential surfaces are labeled according to the diffraction 
channels they asymptotically connect to. The detuning is <5 = 
25d. 



are either repulsive or attractive. For each type, the po- 
tentials form an asymptotically equidistant series with 
a separation equal to twice the Doppler shift. [Strictly 
speaking, the levels are not equidistant due to the 're- 
coil term' quadratic in Q, but this is typically a small 
correction at grazing incidence.] For a positive detun- 
ing, repulsive potentials connect to the ground state 
channels discussed above for the 'one-level atom', while 
attractive potentials asymptotically connect to excited 
state channels {i' odd). Inside the evanescent wave, we 
observe that repulsive and attractive potential surfaces 
approach each other, although without crossing. This is 
because the ground and excited states are coupled by 
multiphoton transitions, so-called Doppleron resonances 
48 . These 'avoided crossings' play a key role in the 
diffraction mechanism, as is discussed now. 



3.2 Coupled wave analysis 

To solve the Schrodinger equation with the Hamiltonian 
(pO|), it is convenient to use a mixed position-momentum 
representation. This takes advantage of the discrete mo- 
mentum change along the x-direction, while the dynam- 
ics normal to the grating is treated in the position rep- 
resentation. The wavefunctions are thus expanded in 
Fourier series with respect to the x-coordinate: 



V even 

ipe{x,z) = ^ ilj^{z)expikxuX 



(23) 
(24) 



Ddd 



In this representation, the coupling term in the Hamil- 
tonian (^o|) becomes an infinite-dimensional matrix V,^^ , 
and the stationary Schrodinger equation takes the fol- 
lowing form 

'^V'' + E^''mW^. = ?#^. (25) 



2M dz^ 



2M 



A direct numerical solution of these equations is in prin- 
ciple possible and has been performed Re- 
cently, also some analytical work has been reported on 
the reflection problem |Q . 

3.2.1 Adiabatic potentials. In order to simplify 
the solution and to gain qualitative understanding, 
Deutschmann et al. diagonalize the coupling matrix 
Vufi{z) in the coupled- wave equations (p5|). Its position- 
dependent eigenvalues Wu(z) are called the 'adiabatic' 
or 'dressed state' potentials, where the latter term 
draws on the analogy to a two-level atom in a standing 
laser wave 



45 



An example of the adiabatic potential 
surfaces is given in Figj^. 

One observes in this figure that far from the evanes- 
cent wave, the potential surfaces are of two types that 



3.2.2 Diffraction mechanism. Consider a ground state 
atom that enters the stationary evanescent wave in the 
channel i/ = and follows the potential Wo{z). In order 
to have diffraction, it has to make a transition to a dif- 
ferent potential surface, otherwise it would be reflected 
in the v = channel or hit the surface. Such a transi- 
tion occurs if the atom cannot follow adiabatically the 
potential surface Wo{z). This is indeed possible because 
the transformation that diagonalizes the coupling ma- 
trix Vi,f^{z) depends on the position z and does not com- 
mute with the kinetic energy operator in the Schrodinger 
Eq.([25|). The dressed levels thus become nonadiabatically 
coupled. In order to compute the nonadiabatic transi- 
tions between the adiabatic potentials, Deutschmann et 
al. exploit the circumstance that nonadiabatic couplings 
are only important between levels that are 'close' in en- 
ergy. As can be seen from FigJJ, this implies that the 
couplings are spatially localized around avoided cross- 
ings of the adiabatic potentials. In their vicinity, one 
may make a two-state approximation and use a gener- 
alization of the Landau-Zener formula to compute the 
wavefunction amplitudes in the two dressed states after 
the atom has passed the avoided crossing (see for 
details).^ 

We are thus led to the following picture for the 
diffraction process: the incoming atomic wave propa- 
gates through the dressed levels and is split and re- 
combined in the avoided crossings. This creates a num- 
ber of partial waves that propagate either towards the 
dielectric surface or escape into vacuum, after having 
been reflected from a repulsive potential. The station- 
ary evanescent wave is hence physically equivalent to 
an array of beamsplitters (avoided crossings) and mir- 



The picture presented so far is only accurate in a semi- 
classical regime where a generalized Landau-Zener theory for 
the avoided crossings applies. It becomes questionable for 
atoms with a low kinetic energy because of wave-mechanical 
effects. In this regime, Deutschmann et al. numerically solved 
a multi-component Schrodinger equation with appropriate 
boundary conditions. 
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Fig. 5 Typical path of an atom through the adiabatic po- 
tentials that leads to diffraction into the v = —2 ground 
state channel. The four circles mark avoided crossings. The 
optimized detuning is chosen: 5 = 2Sd- 

TOTS (repulsive potentials), and the diffraction pattern is 
determined from the beamsplitters' efficiencies and ttie 
'optical path lengths' between the splitters and mirrors. 
The atoms traverse this array much like balls on a Gal- 
ton board with the difference that all the different paths 
are in principle taken at the same time, their amplitudes 
having to be added to obtain the amplitude for diffrac- 
tion into a given order. 

The performance of the reflection grating is limited 
by the fact that the evanescent wave has to achieve both 
diffraction and reflection of the wavefunction. These 
tasks are actually incompatible: for vanishing modula- 
tion e = (single running evanescent wave), the ground 
state potentials are repulsive, yet there are no avoided 
crossings with attractive potentials. For maximum mod- 
ulation e = 1 (a pure standing wave) on the other 
hand, strong avoided crossings lead to adiabatic poten- 
tials widely separated in energy and essentially flat; the 
incoming atoms traverse the evanescent field and hit the 
surface. The most efficient operation of the diffraction 
grating is obtained for an intermediate contrast. Another 
important condition is to achieve an optical potential 
comparable to the Doppler shift: this assures that the 
potential surfaces for the ground and excited states ac- 
tually cross. The available laser power being limited, one 
has to reduce the Doppler shift Qvxi to meet this condi- 
tion [^H]. At optimized operation, the diffraction proba- 
bilities are still small, however: recall that in an avoided 
crossing, the atom makes a transition to an adiabatic 
potential of excited state character. In order to end up 
finally in a ground state diffraction channel, it must pass 
several crossings. Deutschmann et al. argue that at least 
four crossings are involved, as shown schematically in 
Fig.|. 

Since one of these beamsplitters (marked B) is passed 
twice, first in refiection and second in transmission, it 
is optimized with a 50 : 50 splitting ratio. Assuming 
a perfect reflection at the innermost crossing (C) and 
a similar splitting ratio for the two outer crossings (A, 
D), one is led to a diffraction probability of the order 




detuning {S/W^ S^) 

Fig. 6 (courtesy of R. Deutschmann, Ref. Jl^,^) Predic- 
tions of the two-level model for the populations as a func- 
tion of the laser detuning (in units of the 'recoil frequency' 
5k = fiK? /2M). Thin solid curve: specular {v = 0) order; 
thick solid curve: even nonzero orders; dashed curce: odd or- 
ders. The Doppler shift is taken equal to Sd ~ 6.52 x 10^ 5k- 
For a detuning S = 25d, the even diffraction orders u — ±2 
are maximized to a total population of about 6 %. The 
Stueckelberg oscillations are averaged over by convolution 
with a finite energy width for the incoming atoms. 

of (1/2)^ ~ 6 % for the v = —2 ground state channel. 
This value is indeed typical for the results obtained in 
numerical calculations as shown in Figj6[ Diffrac- 
tion being optimized, the specular population {ly — 
channel) is only around 20 % because in a number of 
crossings, the atoms are transmitted to the prism surface 
and lost. The predicted diffraction probabilities also con- 
tain delicate oscillating features, so-called Stueckelberg 
oscillations, that are due to interferences between differ- 
ent paths through the potential surfaces (not shown in 
Fig.|). 

3.3 Comparison to one-level model 

The picture developed by Deutschmann et al. seems so 
largely different from the one-level atom that a compar- 
ison of the two models is desirable. 

First, why do these models make different predictions 
for atomic diffraction at grazing incidence? We have 
seen that in the two-level model, the diffraction mech- 
anism crucially depends on a resonant coupling to the 
excited state playing the role of an intermediate state 
for diffraction into a ground state channel. It is obvious 
that the one-level model cannot reproduce this resonant 
mechanism, since it eliminates the excited state right 
from the start, assuming it to be off-resonant. 

While this answers the first question, it leads us to a 
second one concerning normal incidence: how can we in- 
terpret the efficiency of diffraction at normal incidence 
in the adiabatic potential picture? At normal incidence, 
the Doppler shift is very small, and the adiabatic poten- 
tials are hence quite close (with a separation of the order 
of the recoil energy). A simple estimation shows that 
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for typical atomic velocities, there are strong nonadia- 
batic couplings between the dressed states, even outside 
avoided crossings [i^ . The atomic wavefunction is cou- 
pled to other diffraction channels all over the evanescent 
wave, and the idea of spatially localized beamsplitters, 
that is used in the scheme set up by Deutschmann et al., 
fails in this regime. By contrast the one-level model does 
take into account such a spatially extended coupling. In 
the thin phase grating approximation, e.g., this is done 
by accumulating the phase due to the modulated evanes- 
cent wave potential (that provides the coupling) along 
the path the atom follows during reflection. 

Note finally that the preceding argument also allows 
one to interpret why the diffraction efficiency is cut off 
at grazing incidence in the scalar model: with increasing 
Doppler shift, the adiabatic potentials for ground-state 
diffraction channels become more and more separated, 
and the probability of nonadiabatic transitions decreases 
exponentially with the level separation [Kazantsev et al. 
discussed a similar result in 1980 for two- level atoms 
and a transmission grating ||4^]. The atom hence stays 
adiabatically in its initial repulsive dressed level, and no 
transitions to other (ground state) diffraction channels 
take place. 

3.4 Comparison to experiment 

We have to distinguish between two regimes depending 
on the value of the detuning compared to the Doppler 
shift. 

The two-level model predicts optimized diffraction 
for a detuning 5 = 25d (see Fig.^) . Ground and ex- 
cited states are then efficiently coupled by a low-order 
Doppleron resonance. Recall that in this resonance, the 
atom absorbs a number ^ -I- 1 of photons from one run- 
ning evanescent wave and emits / photons into the other, 
counterpropagating, wave. The resonance condition for 
this process depends on the Doppler shift and reads 
5 = {21 + 1)5d in weak fields |48j; it is shifted towards 
smaller detunings in the presence of light-shift potentials 
0. Optimum diffraction is mediated by a three-photon 
Doppleron (the avoided crossing between the channels 
V — 0, 3 in Fig.^ whose 'effective Rabi frequency' is 
comparable to the optical potential. Once the atom has 
been excited, however, it is subject to an attractive op- 
tical potential and likely to hit the dielectric surface 
and be lost from the incident beam. As discussed above, 
this is the dominant process for a pure standing wave 
where the adiabatic potentials are widely split at avoided 
crossings. The two-level model is thus able to describe 
the Doppleron resonance experiments of the Paris-Nord 
(Villetaneuse) group ji^ and of the Canberra group 
where the evanescent wave's reflectivity was measured 
as a function of detuning. 

If the detuning is much larger than the Doppler shift, 
the two-level model predicts a vanishing diffraction efS- 
ciency (see FigJ|). This regime corresponds however to 



the diffraction experiments in Bonn [gj and in Paris- 
Nord (Villetaneuse) There, an atomic beam was split 
upon reflection from a partially stationary evanescent 
wave. Both groups checked that the angle of the sec- 
ondary beam was in accordance with the kinematics of 
a diffraction process. While the Bonn group observed an 
efficiency of a few percent, in apparent agreement with 
the two-level theory, the Villetaneuse group reported a 
splitting of up to 60 : 40 for the v = 0, —2 channels. It 
is impossible to explain these experiments with the two- 
level theory because for a detuning 5 much larger than 
the Doppler shift, the Doppleron resonances yield a neg- 
ligible diffraction efficiency. This was observed numeri- 
cally by Gordon and Savage Il9| and by Deutschmann 



21 , and we would like to give here an analytical esti- 
mate for the Doppleron coupling. 

For the experimental parameters of the diffraction 
experiments [^^, the relevant Dopplerons involve a 
large number of photons {21 -\- 1 between 10 and 20). 
In this regime, eliminating the 21 intermediate states, 
one finds an 'effective Rabi frequency' of order M,Bn, 



1 



fdEe 



V hS 



21-1 



« Vo{zc) (26) 



where Zc is the position of the avoided crossing. This 
Rabi frequency is negligible compared to the optical po- 
tential Vo {zc) for a detuning S larger than the one-photon 
Rabi frequency dEe~'^^'' /Ti. In the adiabatic potential 
picture, the avoided crossing becomes a nearly perfect 
crossing, and the atom follows its repulsive potential as 
if the attractive potential did not exist. The 'beamsplit- 
ter' is thus perfectly 'transparent', no splitting and hence 
no diffraction occur. 



3.5 Summary 

The diffraction of two-level atoms at grazing incidence 
is possible if the atom makes a sequence of transitions 
that involve, at an intermediate stage, dressed states 
of excited character. This mechanism is effective at low 
detunings where Doppleron resonances provide an effi- 
cient coupling. For a pure standing evanescent wave, it 
leads to a reduced number of reflected atoms because the 
excited state is transmitted down to the dielectric sur- 
face. At the large detunings typical for grazing incidence 
diffraction experiments, however, Doppleron resonances 
become negligible, and these experiments call for another 
model. This has been realized by the Bonn group who 
published, shortly after the two-level theory, a proposal 
for an atomic beamsplitter that involves transitions be- 
tween ground state Zeeman sublevels [^,^. The Can- 
berra theory group noticed that Zeeman sublevels are 
necessary for evanescent wave diffraction when they ob- 
served that a numerical integration of the coupled-wave 
equations (^) gave negligible diffraction for the param- 
eters used in the experiment |16|,[f9|]. 
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4 Model with multiple ground state sublevels 

We now review the diffraction theories that take into ac- 
count the Zeeman (magnetic) degeneracy of the atomic 
ground state. We first outhne the corresponding mecha- 
nism and then summarize the theoretical work done so 
fax. 

4-.1 General 

Atomic diffraction occurs when the atom absorbs a pho- 
ton from one running evanescent wave and emits one into 
the other counterpropagating wave. This 'Raman tran- 
sition' corresponds to a momentum transfer and leads 
to the splitting of the atomic beam in momentum space. 
At grazing incidence, the running evanescent waves that 
form the diffraction grating acquire a frequency differ- 
ence in the frame of the atomic beam that equals twice 
the Doppler shift. If the Raman coupling connects the 
same ground state sublevel, it is hence far off resonance 
and no efficient population transfer may take place. But 
the Raman transition may also connect different Zeeman 
sublevels and become resonant if the sublevel degener- 
acy is lifted. This happens, e.g., in a static magnetic field 
(Zeeman effect) or in a suitably polarized light field. The 
first scenario has been explored for a running evanescent 
wave by Deutschmann et al. , in view of building an 
atomic beamsplitter. In the context of the evanescent 
wave diffraction grating, the second scenario was stud- 
ied numerically in the Canberra theory group [p^ and 
analytically in the Orsay group j20|. 

4-2 Physical picture 

For a basic understanding of the diffraction mechanism, 
let us consider the limit of low saturation and a detuning 
large compared to both the Doppler shift and the nat- 
ural linewidth. The excited state manifold may then be 
eliminated adiabatically, and for a ground state of angu- 
lar momentum Jg, the atomic wavefunction is described 
by the 2Jg + 1 components ipm{x, z), m = — Jg, . . . + Jg. 
It is subject to an optical potential V{z) whose matrix 
elements are of the form [we suppose that Doppler shift 
6d and Zeeman shift are negligible compared to the de- 
tuning S] 

q.q' ,me 

X ( Jg, m; 1, q\ Je,me) ( Je, me\ Jg, m'; 1, q') (27) 

where a product of Clebsch-Gordan-coefhcients appears 
on the rhs, and the electric field is expanded in the usual 
spherical basis with coefficients Eg, q = — 1,0, +1. The 
optical potential couples different Zeeman sublevels if 
the field is not in a pure polarization state with respect to 



this basis. The matrix ( |27| ) then contains both diagonal 
and nondiagonal elements. The diagonal elements lead to 
the light-shift of the Zeeman sublevels, that are coupled 
by the nondiagonal elements. 

One may now proceed along similar lines as for the 
two-level theory, and compute the adiabatic potentials 
for the diffraction channels. These channels are labeled 
by the diffraction order 1^ = 0, ±2, . . . and the magnetic 
sublevel m. Without a magnetic field and for fixed ly, 
the magnetic sublevels are asymptotically degenerate, 
while different diffraction orders v are separated by twice 
the Doppler shift (at grazing incidence). An example is 
shown in Fig.|^(a) for a Jg — 1/2 ground state and a 
single running evanescent wave with TM polarization 
(magnetic field vector perpendicular to the optical plane 
of incidence). To a quite good approximation, this wave 
is in fact cr" circularly polarized provided the laser beam 
is incident far beyond the critical angle. [The quantiza- 
tion axis is parallel to the magnetic field vector.] The 
degeneracy of the sublevels is lifted by the light field be- 
cause the diagonal elements of the optical potential ma- 
trix ( |27| ) differ in magnitude for the sublevels m = ±1/2 




V 



Fig. 7 (a) Potential energy surfaces for a J = 1/2 — > Je = 
3/2 transition in a single (nearly) circularly polarized evanes- 
cent wave. The potentials cross where the difference between 
the light-shifts compensates for the Doppler shift. At this 
position, a linearly polarized copropagating wave leads to a 
resonant Raman coupling to the v = —2 diffraction channel 
(see inset). The kinetic energy ^mv^^ of the incident atoms 
is indicated by the thin horizontal line. 

(b) Schematic representation of the diffraction process in 
terms of a Michelson interferometer. The level crossing of 
(a) is represented as a beamsplitter, and the turning points 
in the repulsive potentials as mirrors. 
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[see inset of FigJ^(a)]. As a consequence, the potentials 
for the diffraction orders v — 0,-2 cross when the dif- 
ference between the light shifts (— 1/2|V^(2;)| — 1/2) and 
(+l/2|y(z)|+l/2) is larger than twice the Dopplcr shift. 
For a pure TM polarization, this is an exact crossing 
because there is no Raman coupling (the optical po- 
tential (^7|) has no off-diagonal elements). A coupling 
can be provided if one adds a second evanescent wave 
with linear (Ti?) polarization: starting from the sublevel 
m = +1/2, the atom absorbs a a~ polarized photon from 
the strong counterpropagating {TM) wave and emits a 
stimulated photon with tt polarization into the weak co- 
propagating {TE) wave. The atom thus ends up in the 
m = —1/2 substate of the ly = —2 diffraction channel 
[Fig.0(a) , inset] . In the presence of the second evanescent 
wave, the adiabatic potentials form an avoided crossing 
at the position of the circle in Fig|^(a). There, an in- 
coming wavef unction in the — 0,m = +1/2} channel 
is split in two parts that are subsequently reflected from 
their respective repulsive potentials and recombined af- 
ter the second passage at the crossing. Note that for this 
particular model, the evanescent wave realizes a 'Michel- 
son interferometer' with a single beamsplitter and two 
mirrors, as shown schematically in Fig.0(b). 

4-. 3 Predictions 

4-3.1 Numerical calculation. In Ref . [p^ , the Canberra 
theory group reported a numerical solution of the 
time-dependent Schrodinger equation for an atom with 
Zeeman-degenerate ground and excited states with Jg = 
Jf. = 2. They studied a situation close to the exper- 
iment a diffraction grating formed by two evanes- 
cent waves and an unpolarized atomic beam (all mag- 
netic sublevels equally populated). The population in 
the V — ^2 diffraction channel is found of the order of 
15 % if the optical polarization vectors are tilted (dif- 
ferent combinations of TE and TM polarizations), in 
qualitative agreement with the experiment. Similar to 
the magnetic beam splitter [|2^, only a single Zeeman 
substate is involved in the multilevel diffraction mech- 
anism. This suggests that a transfer efficiency close to 
100 % could be possible with spin-polarized atoms. 

4-. 3. 2 Landau-Zener theory. In the following, we con- 
centrate on analytical results obtained for the Jg = 1/2 
ground state discussed above. The physical picture ex- 
plained for this model can be translated into a simple 
theory when the avoided crossing is treated by means 
of the Landau-Zener model for nonadiabatic transitions 
Assuming that the atom moves through the cross- 
ing with a constant velocity (fixed by energy conserva- 
tion), the Landau-Zener formula allows one to compute 
the probability amplitudes for the two potentials after 
the crossing. One thus obtains the reflection and trans- 
mission amplitudes for the beam splitter in Fig.||(b). 
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Fig. 8 (courtesy of C. Henkel, Ref. po|] ) Predictions of the 
Landau-Zener model (solid line) and the distorted-wave Born 
approximation (points) for the diffraction of a Jg — 1/2 atom 
from an evanescent wave diffraction grating with a weak po- 
larization gradient. The diffracted population is shown as a 
function of the incident velocity component v^i (in units of 
the 'recoil velocity' Tin/M). The intensity \Ete{z = 0)^ of 
the TE polarized wave is equal to 6 x 10"* the TM intensity. 



Calculating the phase shifts between the beamsplitter 
and the turning points in the WKB approximation, we 
end up with the populations of the diffraction channels 
|i/ = 0,TO = +1/2) and = -2,m = -1/2). 

The result of such a calculation is shown in FigJ| 
(solid line), as a function of the incident velocity Vzi 
perpendicular to the grating. The incoming atom is in 
the substate m = +1/2, the diffracted atom being in 
m = —1/2. The diffracted population shows a maximum, 
and this happens if the position of the avoided cross- 
ing coincides with the atomic turning point. Since the 
atom then spends a long time in a region of resonant Ra- 
man coupling, this could have been expected.0For larger 
energies, the atomic wave explores both potential sur- 
faces [cL Fig.^(a)], and the diffraction populations show 
Stueckelberg oscillations due to the interference between 
the two paths. The diffracted population eventually de- 
creases because the atomic velocity becomes larger and 
larger at the avoided crossing, reducing the splitting effi- 
ciency. The coincidence of the atomic turning point with 
the avoided crossing allows us to formulate the following 
condition for optimum diffraction: 



{-l/2\V{z,)\-l/2) E,, + 2h6D 



-l/2|y(zc)| + l/2) 



E. 



(28) 



Due to the common exponential variation of the optical 
potentials, the Ihs is in fact independent of Zc and only 
depends on the polarization of the evanescent field. Re- 
calling that the atom must be reflected from both optical 
potentials, we recover the prediction of Dcutschmann et 

^ If the atomic velocity is precisely zero at the crossing, 
Landau-Zener predicts in fact that the atom adiabatically 
follows its initial potential surface, leading to zero diffrac- 
tion. A very small deviation from this condition is sufficient, 
however, to produce a large transition probability and hence 
efficient diffraction. 
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al. jl^ tliat a liglit-shift comparable to the Doppler shift 
is necessary for efficient diffraction. 

4-.3.3 Distorted-wave Born approximation. Although 
the Landau-Zener model indicates that diffraction is op- 
timized if the atomic turning point coincides with the 
avoided crossing, this circumstance invalidates the model 
itself, because in the vicinity of the turning point one 
cannot describe the atom by a classical particle with a 
given velocity. We may however use an alternative ap- 
proach based on the distorted-wave Born approximation. 
The Raman coupling is then treated as a perturbation 
that induces a coupling between the wavefunctions in 
the initial and final potentials. Note that these wavefunc- 
tions have asympotically different kinetic energies (in the 
normal direction) due to the Doppler effect. Analyzing 
the matrix element in Fermi's Golden Rule in analogy 
to the Franck-Condon factors familiar from molecular 
physics, we expect diffraction to be most efficient when 
the classical turning points of both wavefunctions coin- 
cide, since their amplitude is maximum there.r It is easy 
to see that this condition is in fact identical to (28|). Fig.|| 
shows the result of the distorted-wave Born approxima- 
tion (points). An optimum diffraction is indeed obtained 
close to the (unphysical) maximum of the Landau-Zener 
model (solid line), and both models are in good agree- 
ment for incident energies above the optimum. 



4-.4 Comparison to experiment 

The numerical calculation of the Canberra the- 
ory group is in fair agreement with the Bonn exper- 
iment 1^, given the number of simplifications in the 
theory (evanescent wave of infinite size, omission of 
van der Waals interactions with the dielectric surface). 
For a more detailed comparison, experiments with spin- 
polarized atomic beams and well-controlled optical po- 
larizations would be very useful. The theory indeed pre- 
dicts a quite efficient population transfer (100 % do not 
seem excluded in principle), with promising applications 
for atomic interferometry. 

We also note that experiments at large angle inci- 
dence may be 'simulated' in a normal incidence geom- 
etry, by introducing a frequency difference between the 
two counterpropagating evanescent waves, in a way sim- 
ilar to early diffraction experiments with bichromatic 
standing evanescent waves 1^,0] . Experiments in this di- 
rection have been performed in the Orsay group and 
confirm the multilevel diffraction theory outlined above 
pl[. For a quantitative comparison, however, the van 



This picture yields still another interpretation of the graz- 
ing incidence cutoff in the one-level model: since the initial 
and final wavefunctions are subject to the same light shift, 
their turning points can only coincide if they have the same 
energy, which is impossible due to the Doppler effect. The 
Franck-Condon overlap is thus far from its maximum value. 



der Waals interaction and, possibly, losses from sponta- 
neous emission have to be taken into account. A detailed 
discussion will be published elsewhere. 



5 Conclusion 

The diffraction of neutral atoms from a stationary 
evanescent wave has remained, since its proposal in 1989, 
a fascinating challenge, both experimentally and theo- 
retically. In the last decade, important breakthroughs 
have been achieved: detailed theoretical predictions for 
both grazing and normal incidence as well as successful 
experimental observations. Despite the conceptual sim- 
plicity of the setup, the diffraction mechanism turned 
out to be quite subtle. In close contact with the experi- 
mental efforts, theory has now evolved towards a picture 
where Raman transitions between magnetic sublevels of 
the atomic ground state play a crucial role. This fea- 
ture endows atomic diffraction with a richer structure 
than, e.g., light diffraction |5|], and it also allows one 
to construct tunable and nearly lossless atom-optical 
beamsplitters, using suitably polarized light fields far 
off resonance and/or magnetic fields. The field is still 
active and one observes an increased interest in specifi- 
cally tailored light fields, using microfabricated surfaces 
42,4^, for grazing incidence reflection gratings. At nor- 
mal incidence, it has become apparent that atoms mov- 
ing in an evanescent wave provide a sensitive probe of 
van der Waals like surface interactions The combi- 
nation with multilevel diffraction mechanisms suggests 
to probe this interaction with interferometric resolution 
and for well-defined Zeeman sublevels. 

Another direction of current research may be termed 
'coherent atom optics' where one studies the motion of a 
high-density (and ultimately Bose-Einstein condensed) 
atomic gas in laser fields. For instance, the evanescent 
wave mirror combined with gravity could be used to 
build an atomic resonator 



55 1 where diffraction may 



serve as a convenient output coupler. At high densities in 
this resonator, one has to take into account the large re- 
fraction index of the trapped gas, modifying the evanes- 
cent wave . A second example are atomic waveguides 
in hollow fibers: 'coating' the walls with an evanescent 
wave, efficient guiding has already been demonstrated 



157 



0^,^. Grazing incidence diffraction taught us that 
Raman couplings may mix the ground state sublevels 
and transfer large amounts of kinetic energy perpendic- 
ular to the wall. This problem will have to be faced on 
the route rowards single-mode atomic waveguides. 
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